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ABSTRACT 

A strong emission line at 2.8935/im discovered by Rubin et al. (2001) in an ISO 
SWS02 spectrum of the Orion Nebula is identified as the 4p 3 P — 4s 3 S° multiplet 
of O I. Line formation is due to de-excitation cascades following UV-pumping of high 
3 S° and 3 D° terms and occurs in the O I zone immediately behind the Hydrogen 
ionization front. This cascade mechanism also accounts for permitted O I triplet lines 
in the Nebula's optical spectrum (Grandi 1975). An escape probability treatment of 
the Oi cascades accounts for the strength of the A2.89^m line and suggests interesting 
diagnostic possibilities for the optical lines. 
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1 INTRODUCTION 

Rubin et al.(2001) have recently reported the detection of 
an unidentified (uid) strong emission line in an ISO SWS02 
spectrum of the Orion Nebula. They also report indications 
of the line's presence on two spectra from the ISO archive 
taken at other locations in the Nebula and strong confirma- 
' tion of the line's reality in a UKIRT long-slit spectrum at 
the first location. Assuming the line arises in the main ion- 
ization zone, they determine the line's vacuum wavelength 
to be 2.89350/im. 

Unidentified weak lines in astronomical spectra are of 
course common. But this line is strong, being a factor of only 
3.6 weaker than the nearby Hi 11-5 line at A2.87/im. Given 
the efforts of many astronomers since the 1930's in studying 
the formation of emission lines in photoionized nebulae, ei- 
ther as diagnostics or as coolants, the discovery of a strong 
uid line at this late date is surprising and potentially of ma- 
jor significance. 

Rubin et al. (2001) report their own attempts to identify 
this line. They note satisfactory wavelength agreement for 
forbidden lines from two ionic species: Cr ill and Fe V, but 
argue against either being the correct identification. In this 
paper, an independent effort to identify this line is reported, 
which was started after seeing their IAU Abstract (Rubin 
et al. 2000). The search for candidates was based on the 
Kurucz-Bell (1995) compilation of lines and gf- values. 

2 REQUIRED EMISSIVITY 

If, following Rubin et al (2001), we assume that the uid line 
originates in the main ionization zone, then its emissivity can 



be estimated from its strength relative to that of the nearby 
H line. Adopting T e = 9000K and N e = 4000cm" 3 as typical 
for the Orion Hn region (Osterbrock, Tran & Veilleux 1992) 
and setting N(H + ) = N e , we find from Storey & Hummer 
(1995) that the Case B recombination emissivity of the 11-5 
Hi line at A2.87^i is 8.8 x 10~ 21 ergs _1 cm~ 3 . Accordingly, if 
the unknown emitting species has the same distribution as 
H + , the uid line's emissivity is 

47r.72.89 M ^ 2.4 x 10~ 21 erg s _1 cm~ 3 . (1) 

Any proposed identification for which the emission origi- 
nates within the H II region must be capable of emitting at 
this rate. 



3 REJECTED CANDIDATES 

In this section, candidates found on the basis of wavelength 
coincidence are reported and reasons given for their rejec- 
tion. Because this work partly duplicates that of Rubin et 
al. (2001), the discussion is abbreviated. 

An excellent candidate on the basis of wavelength is a 
forbidden line of Chromium at 2.89349/im. This was dis- 
carded immediately because of the low cosmic abundance 
of this element. Rubin et al. were more persistent but even- 
tually also concluded against this identification, with the 
stronger argument that three other [Cr ill] lines in the IR 
are not seen in their spectra. 

Another similarly acceptable candidate is a forbidden 
line of Fe v at A2.89343/im. Rubin et al. reject this candi- 
date since the ionization potential of Fe IV (54.8eV) slightly 
exceeds that of He II (54.4eV), implying a negligible rate of 
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creation of Fev ions by photoionization. But, given that Ru- 
bin et al. (2000) had evidently failed to find an identification 
from among the ions expected to be abundant in the H II 
region, the Fe v possibility was pursued further. 

Statistical equilibrium calculations to predict emissiv- 
ities for the rich spectrum of [Fe v] lines were carried 
out using the A- values of Garstang (1957) and the colli- 
sion strengths of Berrington (2001). For T e = 9000K and 
N e — 4000cm -3 , the emissivity per Fe v ion of the transi- 
tion 3d 4 3 G 5 -3d 4 3 F2 4 at A2.89343/im is 0.93 xlO" 20 ergs" 1 
cm"' . Accordingly, to achieve the target emissivity given in 
equation (1), the number density of Fev ions must be ~ 
0.26cm" 3 . But at iV e = 4000cm" 3 and with an Fe abun- 
dance n(Fe)/n(H) = -4.49dex (Seaton et al. 1994), the 
number density of Fe atoms is only ~ 0.13cm" 3 . If further 
evidence against this identification were needed, these calcu- 
lations reveal that there are several other [Fe v] lines in the 
IR with greater emissivities, including two from the same 
upper level as this rejected candidate. 



4 PROPOSED CANDIDATE 

The suggested identification is the O I multiplet 4p 3 P — 
4s 3 S°. The permitted components of this multiplet are at 
vacuum wavelengths 2.89330, 2.89352 and 2.89359/im with 
gf-values of 0.50, 2.51 and 1.50, respectively. 

Two mechanisms of line formation will be considered in 
this section in an attempt to justify this identification. The 
aim is to demonstrate that emission in this transition can 
plausibly account for the maximum observed line brightness. 
No attempt is made to give a definitive treatment of line 
formation. 

4.1 Atomic data 

The atomic model is restricted to the O I triplet system. 
Specifically, all terms in the 3 S°, 3 P, 3 D° and 3 F series are 
included up to principal quantum numbers n = 11, 10, 11 
and 10, respectively. Energy levels are the experimental val- 
ues from the NIST Atomic Spectra Database if available. 
For higher n, we take 



where E x is the ionization potential, Ro is the Rydberg 
constant for Oxygen, and a is a constant for a given series 
obtained by fitting this formula to the highest level reported 
in the NIST database. The resulting atomic model has 84 
levels. Note that sublevels belonging to the same term are 
not consolidated into a single level. 

Oscillator strengths for all the permitted multiplets in 
the above O I model were computed by Butler & Zeippen 
(1991) as part of the Opacity Project and are available in 
the TOPbase archive. Values for the individual transitions 
were derived assuming LS-coupling. 

Radiative recombination coefficients for the individual 
triplet terms for a wide range of electron temperatures have 
been computed by Nahar (1999). Values for the sublevels are 
derived assuming proportionality to their statistical weights. 

Grotrian diagrams for the O I triplet system have been 
published by Grandi (1975) and Przybilla et al. (2000). 



4.2 Recombination cascade 

By analogy with the H emission-line spectra of H n regions 
and planetary nebulae, the most obvious emission mecha- 
nism for this Ol multiplet is via a radiative cascade following 
recombinations of + ions to high triplet levels. 

Because the ionization potentials of Ol (13.618 eV) and 
Hi (13.598 eV) are essentially identical, the emitting volumes 
for O I and H I recombination lines are also identical, namely 
the H II region. Accordingly, on the assumption that it is a 
recombination line, the viability of this identification can be 
tested by comparing the O I multiplet 's emissivity with the 
estimate given in equation (1). 

The emissivity is calculated at N e — 4000cm" 3 and 
T e = 10 4 K, this being a temperature at which Nahar 
(1999) tabulates recombination coefficients. The Oxygen 
abundance is taken to be n(0)/n(H) = — 3.38dex, the B 
star value adopted by Savage & Sembach (1996). 

The calculation proceeds as follows: the rate of recom- 
binations to level i per unit volume is a;n(0 + )n e , where a.i 
is the recombination coefficient for level i. These recombi- 
nations create a cascade down to the ground term. Accord- 
ingly, let Ri denote the rate per unit volume at which level i 
is populated by recombinations directly to i and by cascades 
from higher levels. Then, since the fraction 

Pij = Ajj/y^Au (3) 

i 

of the population of level i decays to level j < i, where Aij is 
the Einstein A- value for the transition i —> j, the emissivity 
of this transition is 

= PijRihvij . (4) 

Moreover, these i — » j decays make the contribution PijRi 
to the quantity Rj . 

The above calculation is started by setting R u — 
a„n(0 + )n e , where u denotes the highest level in the atomic 
model, and then proceeds downwards level by level to the 
ground term, resulting in emissivities for all permitted lines 
in the recombination cascade. 

For the three lines contributing to the 4p 3 P — As 3 S° 
multiplet, the total emissivity is 5.6 x 10~ 24 erg s" 1 cm" 3 , 
a mere 0.23 per cent of the required value. But this is a 
case A calculation in which photons emitted in decays to 
the ground term escape. For case B, which is obtained by 
repeating the above calculation with An = for I — 1 — 3, 
the emissivity increases to 2.3 x 10~ 23 erg s" 1 cm" 3 , still 
only 0.96 per cent of the required value. 

4.3 Pumped cascades 

If emission in the Ap 3 P - 4s 3 S° O I multiplet were possible 
only following recombinations of + ions, then this candi- 
date would now be decisively ruled out. But since the ground 
term of O I is also a triplet, A2.89/im emission can also oc- 
cur via UV pumping of high triplet levels. Accordingly, the 
suggested emission region is the neutral O zone which, be- 
cause of the essentially identical ionization potentials noted 
earlier, is coextensive with the Hi region - see Fig. 1. In this 
region, the far UV but non-ionizing radiation - ie. A > 911 A 
- of 9 1 Ori C penetrates and excites O I atoms from the three 
levels of the ground term to high 3 S° and 3 D° terms, from 
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Figure 1. Geometrical configuration assumed for the radiative 
transfer calculation. Radiation from 9 1 Ori C is incident at angle 
6q on a semi-infinite uniform layer of neutral Hydrogen and Oxy- 
gen. Line emission created in the pumped cascades propagates 
towards the Observer at angle Oi. 



whence radiative decays give rise to emission in the A 2.89 
/xm multiplet. Exactly the same mechanism was invoked by 
Grandi (1975) to explain the permitted O I triplet lines in 
optical spectra of the Orion Nebula. 

For this line-formation mechanism, the H II region is 
not the source of emission; and so the emissivity estimate 
in equation (1) is no longer the relevant target. Instead, 
the proposal must be tested by predicting the line's sur- 
face brightness. In particular, a viable identification must 
account for the maximum brightess measured by Rubin et 
al. (2001), 

/2.89m =0.12 bu , (5) 

where a brightess unit (bu) of 10 -13 erg s _1 cm~ 2 arcsec -2 
is adopted. 

In the remainder of this section, a crude theory is de- 
veloped with the deliberate aim of overestimating the line's 
brightness. If this overestimate were to fall short of this max- 
imum, a more detailed treatment would not be necessary. 

In computing the line's brightness, we assume that the 
UV pumping and subsequent de-excitation cascades occur in 
a semi-infinite uniform slab of neutral O atoms whose sur- 
face z — is coincident with the Hydrogen ionization front. 
The incident pumping radiation from 8 1 Ori C has direc- 
tion cosine /no = cos 8o, and the observer receives radiation 
emerging at fii — cos 9\ . 

Let us consider the pumping line at frequency Vk cor- 
responding to one of the permitted transitions i — > u for 
I = 1 — 3 and u > 4. On the assumption of zero attenua- 
tion in the H II region by either dust or neutral O atoms, 
the intensity at z — of the radiation from 8 1 Ori C is 
Ik — F(v k ), where irF(v k ) is the emergent flux at v k from 
the star's atmosphere. Accordingly, the rate per unit area 
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at which energy in the frequency interval (Av)k centred on 
i>k crosses into the O I zone is 

E k = huo (Ai/) fc (Aw)* . (6) 

Here (Aw)" = iv(R/d) 2 is the solid angle subtended by 9 1 
Ori C from the point considered at the ionization front. Note 
that the contribution of diffuse radiation from the Hn region 
is neglected. 

If we now assume that all the energy in the interval 
(Av) k is absorbed in the slab by neutral O atoms, then 
Ek /hvk is the rate at which the Vk line photons excite level 
u per unit area. Let us further suppose that each of these 
pumping events leads on average to the emission of e%j (<1) 
cascade photons of frequency Vij. It then follows that 

f-OO 

An / j Vij dz = e% V -^E k , (7) 
Jo Uk 

where Amj Vij is the emissivity of the line v^ . 

Because the cascade photons of interest are in the IR 
and are emitted by subordinate lines, we can neglect atten- 
uation as they propagate to the surface of the slab. The 
transfer equation to be solved is therefore simply 

■ J. /o\ 

^~dz~ = 3vii ' ( ) 

where </>„ is the cascade line's normalized emission profile. 
For a uniform O I slab, this profile is independent of z, and 
so the emergent intensity is 



W = 




Now if i? x i? is the element of sky used to define surface 
brightness, then, neglecting extinction between the ioniza- 
tion front and the observer, the line brightness at frequency 
v is I v (ni) i? 2 . Accordingly, if we integrate over the line pro- 
file, the contribution to the brightness of the cascade line 
i — > j from the kth pumping transition is 

, i? 2 f°° 

fii = — / M 3 dz . (10) 
Mi Jo 

Let us now suppose that isotropic macroturbulence with 
characteristic velocity vd is the cause of kinematic broad- 
ening in the O I slab and write (Av)k = WkAvo, so that Wk 
is the dimensionless equivalent bandwidth for absorption by 
the feth pumping line. Then, using equation (7) to eliminate 
the emissivity integral and substituting Ek from equation 
(6), we obtain 

fij = 7^ x — ("j) x e iJ Wk x —F(vk)Av D ■ (11) 
4 /ii a Vk 

Equation (11) gives the contribution to the strength of 
the line i — > j from the cascade driven by the kth pumping 
transition. For the A2.89/im multiplet, these pumping tran- 
sitions comprise all the permitted transitions connecting the 
three states of the ground term 2p 4 3 P to the states belong- 
ing to the terms ns 3 S° with n > 5 and to the terms nd?D° 
with n > 4. Summation of equation (11) over these pumping 
transitions gives the required line strength. 

In the absence of predictions for the equivalent band- 
widths and for the efficiencies with which cascade photons 
are emitted, we make the optimistic estimates that, for all k, 
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Wk = 5 and = 1. In addition, we assume blackbody emis- 
sion by 9 1 Ori C with T e}f = AQQQQK, and take R = II.IRq 
(Lucy 1995). 

Because Rubin et al. (2001) observed maximum line 
brightness at position 1SW, a line-of-sight that passes 9 1 
Ori C with a small impact parameter ol ~ 0.08pc, we set 
d = 0.2pc, which is Wen & O'Dell's (1995) estimate of this 
star's distance above the ionization front. In addition, since 
the line-of-sight with zero impact parameter has /io = Hi, 
we take the orientation factor Ho/Hi — 1- 

Finally, the macroturbulent velocity vd must be esti- 
mated. From their UKIRT spectrum, Rubin et al. (2001) 
measure the resolved FWHM of the A2.89/^m line to be 24 
km s~ . If the line were single, this would give vd = 14.4 km 
s _1 . But allowing for the three unresolved components with 
relative strengths from Section 5, we obtain vd = 12.7 km 
s~\ a value that will be used throughout this investigation. 
Of course, the actual broadening may be dominated by a 
differentially expanding flow, with a much smaller turbulent 
component. Evidence of non-thermal broadening of other 
lines in the spectrum of the Orion Nebula are summarized 
by O'Dell (2001b). 

With these choices of the parameters, the predicted line 
brightness is 



2.89m 



1.48 bu 



(12) 



which is a factor 12.3 greater than the maximum of 0.12 
bu observed by Rubin et al (2001). By thus comfortably ex- 
ceeding observed line strengths, the proposed identification 
survives this initial test. In fact, in the absence of alterna- 
tives, this is already a powerful argument in favour of the 
O I multiplet. 

Note that, according to this crude theory, the bright- 
ness is independent of the Oxygen abundance. This arises 
because of the assumption that all photons within a fixed 
bandwidth are absorbed by the pumping lines. 



5 IMPROVED THEORY 

The simple calculation of Sect. 4.3 is consistent with the 
identification of the A2.89pim line as the Ol multiplet 4p P — 
4s 3 S°. Now an escape probability treatment is developed to 
see if this identification remains viable when the optimistic 
estimates of Wk = 5 and ejj- = 1 are replaced by actual 
calculations. In addition, if the identification is confirmed, a 
more realistic theory will be needed to exploit the measured 
strengths of this and other O I lines formed by the cascade 
mechanism. 



5.1 Incident intensities 

In Sect. 4.3, the incident intensity Ik = B„ k (T e ff), with 
T e f f — 40000K. But the pumping lines are in the wavelength 
range 918-979 A where line blocking by the Lyman series and 
by metal lines is expected to be strong. Accordingly, the first 
improvement is to replace the black body approximation 
by the Kurucz (1979) LTE line-blanketed atmosphere with 
T e ff = 40000K, log g — 4.5 and solar metal abundances. 
This model was previously used to investigate the fluorescent 
excitation of [Nin] and [Fell] lines in the Orion Nebula (Lucy 
1995). 



As a result of this change, the incident intensities for 
lines with A < 930 A are reduced by up to 25 percent while 
those with longer wavelengths are increased by up to 60 
percent. Yet more accurate would be a NLTE model inco- 
porating spherical extension and the stellar wind (Gabler et 
al. 1989). 

5.2 Equivalent bandwidths 

The rate per unit area at which the kih pumping line absorbs 
energy from the incident beam is 



E k 



i v J v dv 



(13) 



where i v is the line's absorption coefficient per unit volume 
and J„ is the beam's mean intensity at depth z. 

On the assumption that the pumping is confined to a 
layer whose thickness is small compared to the distance d of 
8 1 Ori C, J v = I v x (Auj)* /4-7T, where I v , the beam's specific 
intensity, is given by 



h = I k exp[-(k u + . 



+ m v ) — ] 
Mo 



(14) 



Here fc„ is the absorption coefficient of unit slab volume 
due to interstellar grains and m v is the summed absorption 
coefficients of any other pumping lines that overlap the line 
considered, and which therefore compete for photons. 

If we substitute for J„ in equation (13) and integrate, 
we recover equation (6) except that (Av) k is now replaced 
by WkAvo, with equivalent bandwidth w k given by 



Aud / fc„ + l-u + m, 



dv 



(15) 



5.2.1 Dust absorption 

The extinction coefficient for a typical mixture of interstellar 
grains can be derived from standard fits to observational 
data as follows: a line of length s in the O I slab corresponds 
to optical depth t\ = k e x xt s = 0.4£nl0 x A(\), where A(X) 
is the extinction at A in magnitudes. This can be expressed 
in terms of the colour excess E(B — V) by writing A(X) = 
A(X)/A(V) xR v x E(B-V), where R v = A(V)/E(B - V) 
is the ratio of visual to selective extinction. But the colour 
excess is found to be proportional to the H column density. 
Thus, we also have N(H) = n{H)s = -yE(B - V), with 
7 ~ 5.8 x 10 21 atoms cm~ 2 mag -1 (Bohlin, Savage & Drake 
1978). Eliminating s from these formulae, we find that 



kl xt = O.UnlO x ^± x R v x 



A(X) 



n(H) 



A(V) 



(16) 



This is evaluated using the equation A(X)/A(V) = a(x) + 
b(x)/Rv, where x is the reciprocal wavelength in /im _1 and 
the functions a and b are the fits to observational data given 
by Cardelli, Clayton & Mathis (1989). Note that since the 
pumping lines are such that 10.2 < x < 10.9, this appli- 
cation implies a modest extrapolation beyond the observed 
extinction data. The ratio Rv is set = 3.1, the standard 
value for the diffuse interstellar medium (eg, Cardelli et al. 
1989). 

Equation (16) gives the sum of the dust absorption and 
the scattering coefficients. If vj\ is the albedo at A, the ab- 
sorption coefficient is 
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(17) 



In this investigation, we set vj\ — 0.4 for all the far-UV 
pumping lines. This value is suggested by Fig. 8 in Gordon 
et al. (1994), which summarizes the work of these and earlier 
workers. 



5.2.2 Line absorption 

On the assumption that isotropic macroturbulence with 
characteristic velocity vd is the dominant broadening mech- 
anism, the normalized line absorption profile is 



— (- 

Ap d Av d 



Vk f\ 



(18) 



where Avo/vk = vd/c. The resulting line absorption coef- 
ficient per unit volume is 



m e c z 



fin X 



ri{ 

7W) 



x n(0) x 4> v 



(19) 



Here the index £ = 1 — 3 refers to the three levels of the 
ground term, and the index u refers to the pumped levels 
belonging to the series ns 3 S° and nd 3 D°. 



5.2.3 Computed bandwidths 

Before calculating the quantities Wk from equation (15), line 
overlaps must be identified. No overlaps occur for the mul- 
ns 3 S° because the lower levels are well 



tiplets 2p 4 3 L 

separated and the upper levels are single. But for the mul- 
tiplets 2p 4 3 P —> nd 3 D° overlaps do occur because of the 
small fine structure splittings of the upper terms. Specifi- 
cally, the three components 2p 4 3 P2 — > nd 3 £>i,2,3 compete 
for photons, as do the two components 2p 4 3 Pi — > nd 3 D° 2 . 

In contrast to the calculation of Sect. 4.3, which was 
independent of the Oxygen abundance, this improved theory 
depends on n(O) /n(H) through the ratios i v /k v and m v /k v . 
As in Section 4.2, we take n(0)/n(H) = —3.38 dex. 

In addition, the populations of the levels comprising the 
ground term must be specified. Assuming proportionality to 
statistical weights, we set ni/n(0) = (2J + l)/9. In fact, 
these values are in fair accord with the detailed calculations 
of Pequignot (1990) for the expected conditions immediately 
behind the ionization front - i.e., n(H) J> 10 5 cm~ 3 , T ^ 5000 
K, n e /n(H) <, 0.01. 

The equivalent bandwidths Wk obtained from equation 
(15) for all relevant pumping lines are plotted in Fig. 2. 
These widths are all less than the optimistic estimate of 5 
used in Sect. 4. 3. Nevertheless, a substantial number of lines 
have Wk <; 2, thus indicating effective pumping. 

The general trends shown in Fig. 2 are readily un- 
derstood. Each 3 S° term contributes three lines and their 
u>k values fall with increasing n due to declining ^/-values, 
which imply greater losses to dust absorption. For the six 
lines from each 3 D° term, there is the additional effect of line 
overlap. As a result, overlapped lines with small (//-values 
are starved of photons and thus make little contribution to 
the pumping. This effect explains why, for each n, three of 
the n 3 D° components have Wk 1.3. 
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Figure 2. Equivalent bandwidths Wk from equation (15) plotted 
against the excitation energy of the pumped level. Results arc 
shown for all the pumping transtions contributing to emission in 
the Ol A2.89 fim multiplet. Large filled circles identify transitions 
whose upper levels are n 3 S° for n = 5 — 11. Small filled circles 
identify transitions whose upper levels are n 3 Dj 23 forn = 4— 11. 



5.3 Escape probabilities 

In cascade calculations, one of two extreme cases is usually 
considered: case A, in which photons emitted in decays to 
the ground term escape the nebula; or case B, in which all 
such photons are re-absorbed by the emitting transitions. In 
contrast, here we use escape probabilities to treat interme- 
diate circumstances. 

For subordinate transitions - i.e. u — > I for I > 3 , 
the probability of re-absorption is negligible, and so we take 
their escape probabilities [3 u e = 1. But for decays to the 
ground term 2p 4 3 P, the possibility of re-absorption needs 
to be considered. If a photon of frequency v is emitted in the 
decay to level 2p 4 3 Pj, the probability of it being absorbed 
by a dust grain is k v /(k v + i v + m„), where the notation 
is that of Section 5.2. Thus, by averaging over the line's 
emission profile, we find that the escape probability is 



k„+£„+ m v 



dv 



(20) 



for t = 1 - 3 and u > 4. 

If a photon emitted in a decay to the ground term is 
not absorbed by a dust grain, it is here assumed to be re- 
absorbed by the emitting transition. But this neglects two 
effects. First, there is the additional possibility of escaping 
by crossing back into the H II region - i.e. z < in Fig. 
1. Second, there is the possibility - for decays from the 3 D° 
terms - of being absorbed by an overlapping component, thus 
transferring excitation to a different sublevel of the same 3 D° 
term. 

The first possibility is neglected to avoid a z— dependent 
cascade calculation, a complexity not warranted at this 
point. The second possibility is neglected since excitation is 
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Escape probabilities 
due to abs by grains 



1.02 1.04 1.06 1.08 

E„ pp . r (10 s cm >) 

Figure 3. Escape probabilities /3 U £ from equation (20) plotted 
against the excitation energy of the emitting level. Large filled 
circles identify transitions whose upper levels are n 3 S° for n = 
5 — 11. Small filled circles identify transitions whose upper levels 
are n S D° 2 3 for n = 4 — 11. Line absorption coefficients are 
computed on the assumption of Doppler broadening by thermal 
motions at T = 1Q 4 K. 



not lost, merely transferred. As such, this assumption should 
have little effect on predicted line strengths. (In fact, a code 
incorporating this interlocking effect has been written, and 
this statement is confirmed.) 

In computing escape probabilities with equation (20), 
fc„ is calculated exactly as in Section 5.2.1. But the calcu- 
lation of the line absorption coefficients i v and m v and the 
emission profile (f> v departs slightly from Section 5. 2.2 in that 
the macroturbulent velocity vd is replaced by yj2kT /mo, 
the thermal velocity of the Oxygen atoms. This allows for 
the possibility that the length scale on which photon re- 
capture occurs is less than that of the bulk motions that 
determine the broadening of the A2.89^im line. The temper- 
ature immediately behind the ionization front is taken to be 
T = 10000K (Esteban, Peimbert & Torres- Peimbert 1999). 

The resulting escape probabilities are plotted in Fig. 3. 
As for Fig. 2, the general trends are readily understood as 
being due to decreasing (//-values as one ascends the two 
series. But one point of difference is that the D° terms pro- 
vide three points on Fig. 3 as against six on Fig. 2. This is 
because the quantity i v + m v in equation (20) is identical 
for the three overlapping lines nd 3 D°2,3 ~> 2p 4 3 Pz, which 
therefore have identical escape probabilities. Similarly, the 



two overlapping lines nd 3 D° 2 
escape probabilities. 



5.4 Efficiencies 



2p Pi also have identical 



With escape probabilities determined, the efficiencies can 
now be derived from the cascade driven by the kth pumping 
transition. The calculation of the pumped cascades follows 
closely that for recombinations described in Section 4.2. 



Efficiency 
in creating 
2.89 fi photons 



! *. 



1.02 1.04 1.06 1.08 

E t „ mp . a (cm >) 

Figure 4. Efficiencies eij 8g, M from equation (22) plotted against 
the excitation energy of the pumped level. Large filled circles 
identify transitions whose upper levels are n 3 S° for n = 5 — 
11. Small filled circles identify transitions whose upper levels are 
n 3 D° 2 3 for n = 4 - 11. 



The photon absorption rate per unit slab area in the 
feth pumping line is E k /hv k . The quantity E k is derived 
from equation (6) with n — 1 and (Au) k = w k Av D , with 
equivalent bandwidths from Section 5.2. 

If the fcth pumping line corresponds to the transition 
I — » u, then the kih cascade starts at level u with R u — 
E k /hv k . From a level i < u, the fraction of decays that go 
to level j is 



Pij = Aijflij/ ^2 A t e(3ie 



(21) 



where the escape probabilities are from Section 5.3; and 
these decays make the contribution pijRi to Rj. As in Sec- 
tion 4.2, this calculation proceeds downwards level by level 
until the ground term is reached. 

When the cascade terminates, photon emission rates 
PijRij are available for all permitted cascade transitions 
i — > j. But these rates were achieved by absorbing pumping 
photons at the rate R k = E k /hv k . Accordingly, the effi- 
ciency with which the fcth pumping line creates i — » j pho- 
tons is 



e a = PijRi /Rk 



(22) 



Note that because of the cascade's linearity the efficiencies 
efj are independent of fio- 

The efficiencies with which the numerous pumping lines 
create A2.89 am photons are plotted in Fig. 4. The highest 
efficiencies of ~ 0.51 are found for the six components of 
the multiplet 2p 4 3 P — » 5d 3 D°. But most are in the range 
0.1-0.2, a substantial reduction from the optimistic estimate 
°f £ 2.89m — 1 used in Section 4.3. 

Surprisingly low efficiencies are found for the compo- 
nents of the multiplet 2p 4 3 P —> 4d 3 D°, with values in the 
range 0.06-0.13. Inspection of the branching ratios for decay 
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channels from the 4 3 D° levels shows that, despite the small 
escape probabilities for decays to the ground term - see Fig.3 
- these still dominate, accounting for 55-80 per cent of the 
total decay rates from these levels. 

The efficiencies plotted in Fig.3 are not an immediately 
reliable guide to the relative importance of the numerous 
pumping transtions since a line of high efficiency might have 
a low equivalent bandwidth. Accordingly, the contribution of 
each pumped term to the emission of A2.89/im photons has 
also been determined. With their percentage contributions 
in parentheses, the most important terms are: 5 3 D° (29.4), 
5 3 S° (13.7), 6 3 D° (10.6), 6 3 S° (7.7), 4 3 D° (6.9), and 7 3 D° 
(6.4) 

5.5 Predicted intensities 

With the quantities Wk and now calculated, surface 
brightnesses for all multiplets formed in the cascades are 
obtained by summing equation (11) over the pumping tran- 
sitions and over the multiplets' components. The results for 
all multiplets with fx > 0.008 bu are given in Tables 1 and 2. 
Table 1 lists the optical lines, with wavelengths in air given in 
A. Table 2 lists the infrared lines, with vacuum wavelengths 
given in jj,m. The tabulated wavelengths A e // are the mul- 
tiplets' effective wavelengths, computed by weighting each 
component's A by that component's fractional contribution 
to fx. 

The cascades also result in a UV line spectrum extend- 
ing almost to the Lyman limit. These are not tabulated but 
we note that the two strongest UV lines are the 3s 3 S° — 
2p 4 3 P multiplet at A1304.28A and the 4rf 3 D° - 2p 4 3 P mul- 
tiplet at A973.02A. The predicted strengths of these lines 
are 16.7 and 2.11 bu, respectively. 

Of most immediate interest is the predicted strength for 
the 4p 3 P - 4s 3 S° O I multiplet at A2.89^m. From Table 2, 
we find that 

/2.89 M =0.14 bu , (23) 

remarkably close to the maximum value of 0.12 bu measured 
by Rubin et al. (2001) at position 1SW. Since the geomet- 
rical parameters d = 0.2pc and /xo/mi = 1 were chosen to 
represent this point, this outcome is compelling evidence in 
favour of the proposed identification. 



6 DISCUSSION 

In this section, further testing of the proposed identification 
is carried out using the theoretical understanding developed 
in Section 5. In addition, the sensitivity of of the results to 
various input parameters is briefly discussed. 

6.1 Theory-independent estimates 

Because the prediction of about the right strength for the 
Ol multiplet at A2.89/im has required a rather complicated, 
multi-parameter theory, it is of interest that this multiplet's 
strength can be reliably estimated without invoking this the- 
ory. There are two ways of doing this, both of which are 
independent of the line formation mechanism and of the 
distribution of O I along the line-of-sight. 

The first estimate follows from noting that the proposed 



identification has the same upper term (4 P) as the O I 
line at A4368A. Accordingly, whenever the 4368 A line is ob- 
served, we can confidently infer that emission at A2.89/im is 
also occurring and can compute the IR multiplet's strength 
from the formula /2.89m = /o.44/i/l-48 - see Tables 1 and 2. 
This calculation depends only on the A-values and frequen- 
cies of the transitions. 

For a definitive test of the O I identification using this 
accurately known intensity ratio, the measured optical and 
IR line strengths are needed at exactly the same location. 
Unfortunately, this information seems not to be available 
from published spectra. Failing this, the approach adopted 
here is to use the measured strengths of the A4368A line 
at locations near to 1SW to infer strengths for the A2.89/im 
line. This can be done using the data of Esteban et al. (1998) . 

Esteban et al. (1998) report that the A4368A multiplet 
has reddening-corrected brightnesses of 0.036 and 0.060 bu 
at their slit positions 1 and 2. The inferred values for the 
A2.89^m line are therefore 0.024 and 0.040 bu, both of which 
represent averages over the 13.3 "x 2 " slit and refer to po- 
sitions offset from 8 1 Ori C by 45 and 27 ", respectively. 

These inferred values are significantly less than the 
maximum of 0.12 bu observed by Rubin et al (2001) at 1SW, 
a position offset from 6 1 Ori C by 32 ". But this latter mea- 
surement refers to a 0.74 arcsec 2 aperture extracted from 
their long-slit UKIRT spectrum. On the other hand, from 
their ISO S WS02 spectrum with a 14" x 20 " aperture also 
centered on 1SW, they report an average brightness for the 
A2.89/im line of 0.039 bu. This is in satisfactory agreement 
with the spatial averages inferred from the data of Esteban 
et al. (1998) at similar offsets from the UV source. Accord- 
ingly, these theory-independent estimates strongly support 
the identification of the A2.89/wn line as the 4p 3 P — 4s 3 S° 
multiplet of O I. 

The second theory-independent estimate derives from 
noting that other observed optical Ol lines tell us that terms 
higher than 4 3 P also emit. Accordingly, if the cascade from 
such a term results in the emission of the A2.89/xm multiplet, 
this emission's intensity ratio to the observed optical line can 
be computed; and this allows a contribution to the strength 
of the A2.89/im line to be estimated from the strength of 
the optical line. Note, however, that, when the excited up- 
per term is 3 D°, the intensity ratios depend on the relative 
populations of the 3 D° term's sublevels. Also the cascades 
overlap so the contributions to the A2.89/im line cannot sim- 
ply be co-added. Partly because of these complications, this 
means of estimating the strength of the A2.89/im is not pur- 
sued further. But the main reason is that it is more reward- 
ing to use the information in the other optical O I lines to 
test the cascade theory of Section 5. 

These two theory-independent procedures can of course 
be applied to infer the prescence of and to estimate stengths 
for numerous other O I lines in the IR. Since the observed 
optical lines at AA5275 and 5299A imply excitation of levels 
as high as 7 3 D° and 8 3 S°, most of the IR lines in Table 2 
must be present at some level at Positions 1 and 2 of Este- 
ban et al. (1998). But again, rather than proceed with such 
estimates, it is preferable to use the data to test the cas- 
cade theory. If the test is successful, ratios of line strengths 
from Tables 1 and 2 should be reliable and so can be used 
to estimate strengths of IR lines from observed optical lines. 
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Predicted ratios 



Figure 5. Cascade test for permitted O I emission lines in the 
optical spectrum of the Orion Nebula. Open and filled circles 
refer to Positions 1 and 2, respectively, of Esteban et al. (1998) 
and correspond to the indicated wavelength in A. The theoretical 
values arc from Table 1. 

6.2 Cascade test 

As noted earlier, Grandi (1975) has previously demonstrated 
convincingly that pumped cascades account for the permit- 
ted O I triplet lines in optical spectra of the Orion Neb- 
ula. Neverthelees, it is of interest to test the cascade hy- 
pothesis again using the theoretical treatment given here 
together with the modern CCD measurements of Esteban 
et al. (1998). 

The results of the cascade test are shown in Fig. 5 where 
both the observed and predicted line brightnesses are ex- 
pressed as fractions of the summed values for the lines 
AA4368, 5275, 5299, 5513, 5555, 5959 and 6046A. (The line 
A7002A is excluded from this sum because it was observed 
only at Position 2.) Inspection of Fig. 5 reveals a satisfactory 
degree of agreement between theory and observation for all 
but one of the eight observed lines. The exception is the line 
at A7002A. 

A point to note in assessing the success of this test is 
the complexity of this cascade mechanism when compared to 
a simple recombination cascade. In this case, predicted line 
strengths represent the superposition of numerous separate 
pumped cascades whose relative importance is determined 
by dust absorption through its effect on equivalent band- 
widths and escape probabilities. Evidently, there are inter- 
esting diagnostic possibilities to be exploited using accurate 
measurements of the strengths of these and other optical O I 
lines. 

Note that none of the observed optical lines has a 
strength that can be computed via a simple cascade calcu- 
lation from one of the other optical lines. Thus Fig.5 is free 
from this species of degeneracy, which further emphasizes 
its diagnostic importance. 



6.3 Wavelength 

To be accepted, a line identification must be consistent with 
the measured wavelength of the uid line. From their ISO 
SWS02 spectrum, Rubin et al. (2001) infer a rest wavelength 
of 2.89350 ± 0.00003/im, in essentially perfect agreement 
with X eff = 2.893509^m from Table 2. 

However, Rubin et al. derived this wavelength assuming 
that line formation occurs in the main H II region at Vhe.Uo 
— 17 km s~ . But for the proposed identification, line for- 
mation occurs immediately behind the ionization front, for 
which an appropriate velocity is 25.5 km s" 1 (O'Dell 2001a). 
Thus there is in fact a wavelength discrepancy equivalent to 
8-9 km s _1 , significantly larger than the Rubin et al. error 
estimate of 3 km s _1 . But given that a theory- independent 
argument (Section 6.1) places an O I line with about the 
right strength at this slight velocity offset, this discrepancy 
is not grounds for discarding the identification. 

A possible explanation of this velocity offset is that that 
the UV pumping occurs predominantly in neutral blobs en- 
trained in the outflowing ionized gas. 



6.4 Intensity profile 

Rubin et al. (2001) found that the highest surface brightness 
for the uid line is coincident with the brightest part of the 
H n region, and therefore understandably concluded that it 
originates from within the main ionization zone. But the O I 
A2.89^im identification conflicts with this conclusion; and so 
there is the obvious expectation that it will fail to yield a 
satisfactory distribution of surface brightness. 

A crude calculation of the angular distribution of O I 
cascade emission can be made using the 3-D model of the 
Orion Nebula constructed by Wen & O'Dell (1995). Fig.5 
of their paper gives the cross section at position angle 226° , 
which includes the brightest part of the Nebula. To a suf- 
ficient approximation, their profile for the ionization front 
can be modelled as three straight line segments. Taking 500 
pc as the distance of the Nebula, we can use this model to 
compute the distance d in pc from 8 1 Ori C to a point on 
the ionization front. The direction cosines /io and /ii - see 
Fig.l - at this point can also be computed if we add the 
assumption that the ionization front is not tilted in the di- 
rection perpendicular to this cross section. In this way, the 
geometrical factor 



can be derived and used to obtain the intensity profiles of 
the O I cascade lines from the reference line brightnesses 
given in Tables 1 and 2. 

In Fig. 6, the resulting intensity profile for the O I 
A2.89/im multiplet is plotted for position angle 226°. This 
should be compared to Fig. 5a in Wen & O'Dell (1995), 
which plots intensity profiles for [N n], Ha and [O ill] lines 
at this same position angle. This comparison shows that O I 
line has a predicted intensity profile that is qualitatively sim- 
ilar to those of lines emitted by the Hn region. In particular, 
peak intensities at offsets of ~ 30" to the SW are closely co- 
incident, and all profiles show a similar steep decline beyond 
their peaks in the SW. Thus the observed intensity profile 
of the uid line, though suggestive of emission from within 
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Figure 6. Intesity profile for O I A2.89/im emission. This 1-D 
profile is at position angle 226° and corresponds to Fig. 5 in Wen 
& O'Dcll (1995). 

the H II region, is in fact also consistent with line formation 
by pumped cascades in the neutral zone immediately behind 
the ionization front. 

Note that the discontinuities in Fig.6 at offsets of -65 
and +45 "are due to the discontinuities of /Lio//ii at points 
where the linear segments meet. 

6.5 Sensitivity 

In this subsection, the sensitivity of the cascade calculation 
to input parameters is investigated. Given the strong evi- 
dence already presented in support of the O I identification, 
the main interest now is the possible relevance of this line 
and, more importantly, of the optical lines for diagnostic 
investigations. 

6.5.1 Oxygen abundance 

An increase of n{0)/n(H) by O.ldex increases the strength 
of the A2.89pim multiplet by 0.032dex, thus confirming the 
insensitivity to abundance anticipated by the crude theory 
of Section 4. 

With respect to the observed optical lines, the responses 
to this change range from 0.013dex (A5959A) to 0.067dex 
(A7002A). This illustrates a useful differential sensitivity 
that can perhaps be exploited diagnostically. 

6.5.2 Gas-to-dust ratio 

The number density of grains in the O I slab is determined 
by the coefficient 7 in the empirical formula of Bohlin et al. 
(1978) - see Section 5.2.1. Decreasing this number density 
by l.Odex by increasing 7 to 5.8 x 10 22 atoms cm -2 mag -1 
results in a 0.23dex increase in the A2.89/im line. The cor- 
responding changes in the optical lines range from O.lOdex 
(A5959A) to 0.46dcx (A7002A). 

6.5.3 Ratio of visual to selective extinction 

The dust absorption coefficient calculated in Section 5.2.1 
assumes Rv =3.1. This ratio is known to depend on environ- 
ment, with Rv = 5 found in some dense clouds (e.g., Cardelli 
et al. 1989). With this value, the strength of the A2.89Atm 
line increases by 0.044dex. The corresponding changes in 



the optical lines range from 0.018dex (A5959A) to 0.091dex 
(A7002A). 

6.5.4 Microturbulence 

The adopted kinematic model with vd= 12.7 km s _1 has 
been described as macroturbulence because it does not apply 
on the shortest length scales, namely the mean free paths 
of photons emitted in decays to the O I ground term. On 
these scales, thermal broadening has been assumed. If this 
assumption is dropped, and the Doppler parameter in Sec- 
tion 5.3 also set = 12.7 km s _1 , then the kinematic model 
becomes one of microturbulence. 

With this change, the strength of the A2.89/im line 
changes by -0.058dex. The corresponding changes in the 
optical lines range from +0.001dex (A5959A) to -0.30dex 
(A7002A). Notice that this change would seriously worsen 
the already substantial discrepacy for the A7002A line in 
Fig.5. 

6. 5. 5 Temperature of O I slab 

In the calculation of escape probabilities (Section 5.3), ther- 
mal broadening is assumed at T = 10000K. If this temper- 
ature is reduced to 3000K, the A2.89/xm line's strength in- 
creases by 0.035dex. The corresponding changes in the opti- 
cal lines range from 0.005dex (A5959A) to 0.14dex (A7002A). 



7 CONCLUSION 

The aim of this investigation has been to identify the strong 
A2.89^im emission line discovered by Rubin et al. (2001). 
Surprisingly, the identification turns out to be a line from a 
neutral atom (Ol) rather than from an ion with a production 
ionization potential between 13.6 and 54.4 eV as expected 
by Rubin et al. Correspondingly, the line originates not from 
within the H 11 region but from the neutral zone behind the 
ionization front. Because the Orion Nebula is a thin blis- 
ter of ionized gas on the near side of the Orion Molecular 
Cloud (eg, O'Dell 2001b), lines-of-sight through the densest 
part of the H 11 region also intersect the ionization front and 
therefore also the neutral zone where UV pumping of O I 
occurs. This explains why an O I cascade line can have a 
spatial intensity profile that mimics those of emission lines 
originating in the H 11 region (Section 6.4). 

On a memorable previous occasion, puzzling uid nebu- 
lar emission lines also proved to be due to Oxygen (Bowen 
1927), with momentous consequences for astronomical spec- 
troscopy (eg, Osterbrock 1988). In the present case, by con- 
trast, one could argue that the identification is almost totally 
inconsequential given that that the strength of the A2.89/im 
line can be inferred from that of the A4368A line (Section 
6.1). However, the treatment of O I cascades stimulated by 
this uid line reveals interesting diagnostic possibilities, most 
or even all of which can be exploited with measurements 
of the optical O I lines. First, the cascades occur within a 
very narrow zone immediately behind the ionization front. 
For example, with n(H) = 10 5 cm -3 and vd = 12.7km s -1 , 
the length scale corresponding to r = 1 at the line centre 
of the third strongest line pumping the 5 3 D° term is 4 x 
10 _5 pc, equivalent to an angular scale of about 0.02 ". Thus 
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Table 1. Optical O I lines formed in pumped cascades 
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Table 2. Infrared O I lines formed in pumped cascades 
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velocity measurements of the permitted O I lines potentially 
give a sharply defined reference velocity for studies of the 
dynamics of the Orion Nebula. There is also the possibility 
of investigating the dumpiness of O I and therefore H I at 
the point of photoionization by 9 1 Ori C (Section 6.3). In 
view of the abundance of Oxygen and the high ^/-values of 
the pumping transitions, this superb resolution in depth at 
ionization fronts may never be surpassed. 

The second point concerning diagnostic potential is that 
intensity ratios of optical lines are not simply derivable from 
A- values and frequencies (Section 6.2): the astrophysics de- 
termining the relative importance of the numerous indepen- 
dent cascades needs also to be correct; and this leads to 
these lines having usefully different sensitivities to input pa- 
rameters (Section 6.5). 

Yet another possibility is to exploit the geometrical fac- 
tor C - equation (24) - to determine the topography of the 
Orion Nebula's ionization front independently of the proce- 
dure followed by Wen & O'Dell (1995). However, in such an 
investigation, it may be preferable to replace the orientation 
factor /io /mi by an average derived from a statistical model 
of dumpiness. 

These intriguing diagnostic possibilities have not been 
pursued in this paper beyond noting that a substantial evac- 
uation of grains from the cascade zone reduces the A7002A 
line's discrepacy in Fig. 5, and that a microturbulent kine- 
matic model would increase this discrepancy. 
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